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TERMINAL LOCATION 



The present invention relates to a method and an apparatus for locating a g 

• ■ 

mobile terminal in a communications network, particularly but not exclusively 
using multiple sources of information. 



The ability to pinpoint the location of mobile terminals is a desirable feature of 

* 

a mobile telephone network. This is because of the need to provide customer 
services which rely on knowing the whereabouts of users of these services. 
For example, up-to-date local traffic information can be provided to enable a 
user to avoid nearby traffic jams. A user may also wish to know, for example, 
how to get to the nearest pub or restaurant from their present location. 
Clearly the location of the user must be ascertained to within even a few 
metres for this type of service to work. 

Another reason for wishing to know the location of a mobile terminal is so that 
emergency services can locate a caller who is unable to provide an accurate 
personal location themselves. 

It is known in a GSM mobile network to provide the location of a mobile 
telephone in terms of the cell of the network in which the telephone is located. 
Each cell contains one base station and a telephone is only ever in traffic 
communication with one base station at a given time. Hence the location of 
the telephone to an accuracy of the cell area can be <Jetermined simply by 
ascertaining with which base station the telephone is communicating. Such 
methods are known as cell-based location methods. Other methods can be 
combined with the cell identity, such as a trianguiation system, in which-the 
location of a particular mobile phone is calculated using control signals from at 
least the three base stations-closest to it (two of which are located in adjacent 

> 

cells to the cell in which the mobile telephone is located). This system uses 
the assumption that the distance of the phone from a base station is 
proportional to the strength of the signal which the base station receives from 
it, or the time taken for the signal to travel between the phone and the 
respective base station. Thus the position of the phone can be determined by 
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comparing the relative strengths or travel times of received signals between 
the three base stations and thus assessing the distance of the user from each 
base station. The actual location of the user is then obtainable geometrically 
since the location of the base stations is known and fixed. 

In a 3GPP (3 rt Generation Partnership Project) network using a Wideband 
Code Division Multiple Access (W-CDMA) signalling system, it is possible for 
a mobile terminal to be in active communication with more than one base 
station at any one time. This situation is known as "soft handover" and differs 
from (hard) handover in a GSM system, in which a mobile terminal is "handed 
over' from one base station to another as it moves between cells of the 
network. Because of the nature of the soft handover, the above-described 
cell-based mobile location procedures suitable for GSM can not always be 
used in a W-CDMA type signalling system. It would be desirable to provide a 
more reliable way of locating a mobile terminal in this type of signalling 
system. 

In W-CDMA a "softer handover" is defined as well. In the case of "softer 
handover* the antennas of the base stations with which the mobile station is 
communicating are co-located (e.g. they are installed at the same physical 
location). In the remainder of this document, the term "soft handover* will be 
used also to cover the case of "softer handover", and it will be understood by 
those skilled in the art that the invention and the described embodiments 
thereof are applicable to a softer handover situation as well as a soft 
handover situation. 

A problem associated with providing the location of a mobile terminal is that in 
order for the provided location to be meaningful and usable, the accuracy of 
the location provided must also be known. This is because it would be 
pointless, for example, to advise a user of the nearest restaurant to the 
location provided if the user's actual position could be within a range of 
several km. The accuracy can depend on a number of factors such as the 
type of base station antenna {e.g. omni-direction or directional), the -cell size 
{i.e., the extension of the geographical region served by the base station) and 
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the density of network coverage -(number of base stations per square km) in 
the area in which the mobile terminal is located. It would also be desirable to 
know the accuracy of the location of a mobile terminal provided. 

According to the present invention there is provided a method as set out in 
claim 1. 

Further preferred aspects of the invention are set out in the other claims. 

Embodiments of the invention will now be described, by way of example only, 
with reference to the accompanying drawings in which : 

■ 

Figure 1 shows schematically a mobile network divided into service areas 

■ * m 

Figure 2 shows a flow chart for a mobile terminal location procedure in 
accordance with embodiments of the invention 

Figure 3 shows diagrams representing various possible confidence calculation 
methods 

Figure 4 shows the modelling of a cell of the network 

Figure 5 shows two different approaches to evaluation of service area density 

Figure 6 shows geometry for . deriving a circular confidence region and a 
circular arc confidence region according to one embodiment of the invention 

Figure 7 shows a circular confidence region calculation according to a 'Second 
embodiment of the invention. 

Figure 8 shows a polygonal confidence region generated according to a 
variation on the second embodiment of the invention. 
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Figure 9 shows how the location of a mobile terminal is obtained in 
accordance with a third embodiment of the invention 

Figure 10 shows an elliptical and a polygonal confidence region according to 
variations on the third embodiment of the invention 

Figure 1 1 shows a grid for use in another variation on the third embodiment of 
the invention. 

Figure 1 shows schematically part of a Universal Mobile Telephone System 

« 

(UMTS) Public Land Mobile Network (PLMN). The PLMN is indicated by 
reference numeral 1 and it will be appreciated that the network 1 extends 
well beyond the boundary drawn in the figure. Within the network 1 is 
shown a location area (LA) 2. The LA 2 is depicted as being circular but 
this is not necessarily the case in practice. There are in fact a number of 
location areas within the network 1 but the others are not shown. The 
location area 2 comprises a number of service areas 4, in this example four 
(4a^d) as shown in figure 1 distinguished by different shading. Thus a 
Service Area (SA) is defined as a set of cells within a larger Location Area 
(LA). In the example of figure 1 , each SA 4 comprises four cells 6, but an 
SA could comprise a different number of cells. 

■ 

It can be seen that in figure 1 the cells 6 are hexagonal in shape and that 
consequently the shape of each SA 4 is that of four adjacent hexagons. 
This is a well-known approximation to the shape of real network cells. 

The Service Area is identified with a parameter called Service Area 
Identifier (SAI). Thus the SA| can„be_used to identify one or more_ cells 
contained in the same Service Area. Each cell has a unique Cell Identity 
(CI) and has one base station via which mobile terminals can access the 
network 1 . Any cell with which a mobile terminal has an active connection 

■ 

is termed a serving cell. A mobile terminal -can have an active connection 
with more than one cell at any one time. This state is known as "soft 
handover'. 
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The embodiments of the invention which will be described below are 
directed towards calculation of a location estimate and a "confidence 
region* associated with the location estimate of a User Equipment (UE). 
The various embodiments cover the use of single and multiple Cell Identity 
(CI) information. The approaches of the embodiments are suitable for 
application in the context of at least the following two available UMTS 
location methods : 

• Service Area Identifier (SAI) location method 

• Cell Identity and Round Trip Time-(CI+RTT) location method 

The system has a plurality of methods available to it which It can apply to the 
determined location so as to estimate a confidence region in which the mobile 
station could be located to a set probability, 

Before describing the embodiments, these two location methods will be 
explained. 

Service Area Identifier (SAfl Location Method 

In the context of the SAI location method, embodiments of the invention 
perform location calculations based on the fact that the UE is in a certain 
Service Area. This method is intended to calculate the location estimate 
and the confidence region associated with the location estimate of a User 
Equipment (UE) which is reported to be within a known SA of the UMTS 
PLMN of interest. 

The SAI identifying the SA to which a certain «ell belongs is defined by the 
. combination of three codes/identifiers, as follows: 

SAI=PLMN-ld + LAC + SAC 

where: 

PLMN-ld: PLMN Identifier, defined in turn by the combination of 
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- MCC: Mobile Country Code (common to all cells of 
the whole PLMN) 

- MNC: Mobile Network Code (common to ail cells of 
the whole PLMN) 

LAC: Location Area Code 
SAC: Service Area Code 

Any given cell may belong to one or two Service Areas. When it belongs to 
two Service Areas, one is applicable in the broadcast (BC) domain and the 
other is applicable in both the CS (Circuit Switched) and PS (Packet 
Switched) domains of the network. A Service Area in the BC domain 
consists of only one cell, whereas such limitation does not apply in general 
to CS and PS domains. 

The SAl defined in the CS and PS domains (and not the one defined in the 
BC domain) can be used to indicate the location of a UE to the Core 
Network (CN) for LCS (Location Services) purposes. In UMTS, this is done 
by the Serving Radio Network Controller^S-RtMC) which, when requested, 
maps the Cell Identity (CI) of the serving cell into the SAl identifying the 
Service Area to which the serving cell belongs, and sends the SAl to the 
CN by means of the RANAP signalling over the lu interface. 

In general terms, the SAl location calculation algorithm -can be used io 
estimate geographical coordinates and associated confidence region of a 
UE when the only location information available is the identity of cells in the 
Service Area of the serving cell (or cells, when the UE is in^oft handover). 

Cell Identity and Round Trip Time/Ci + R TT) Location Method . 

In the context of a CI+RTT location method, embodiments of the invention 
perform location -calculations based on the fact that when no RTT -or 
RxTxTD measurements <which measurements are in principle available in 
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CI+RTT location method, see below) are used or can be used, information 
from multiple serving cells can be combined to perform location 
calculations. This is possible in a W-CDMA type network due to the soft 
handover functionality mentioned above. 

The CI+RTT location method in UMTS relies on the availability of 
Round Trip Time (RTT) and Rx-Tx Time Difference {RxTxTD) 
measurements. RTT and RxTxTD measurements are introduced in 
UMTS FDD (Frequency Division Duplex) to allow the implementation of 
the CI+RTT location method. 

The RTT is defined as RTT = T Rx UL - Tjx dl , where TDtx dl is the time of 
transmission of the beginning of a downlink dedicated physical channel 
(DPCH) frame to a User Equipment (UE) and Trx ul is the time of 
reception of the beginning (the first detected path, in time) of the cor- 
responding uplink DPCCH (Dedicated Physical Control 
Channel)/DPDCH (Dedicated Physical Data Channel) frame from the UE. 

The RTxTxTD = Ttx 01 -Trx DL is the difference in time between the UE 
uplink DPCCH/DPDCH frame transmission (T Tx UL ) and the first 
detected path (in time) of the downlink DPCH frame from the measured 
radio link (Trx DL ). 

RTTs are measured by the base stations, RxTxTDs are measured by 
the UE. 

By combining a pair of RTT and RxTxTD measurements referred to the 
same base station the distance between the UE and that base station ^can 
be estimated. Such distance estimate is analogous to the distance 
estimate that can be obtained from one Timing Advance (TA) in GSM. In 
this sense, the CI+RTT location method corresponds to the Cell Identity + 
Timing Advance <CI+TA) location method in <3SM. However, -two 
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particular features of UMTS FDD make the CI+RTT method potentially 
more accurate than the CI+TA method in GSM : 

1. The much shorter UMTS chip period as compared to the GSM bit 

t 

period affects the resolution with which a distance estimate can be 
determined from a TA in GSM or from an (RTT, RxTxTD) pair in UMTS. 
One GSM bit period is equivalent to approximately 1100 meters while 
one UMTS chip period is equivalent to approximately 80 meters, thus 
the distance measurements resolution in UMTS is better than in GSM. 

2. In UMTS a UE can be in soft handover. UMTS standards require that 
RTTs and RxTxTDs are measured for each active radio link, thus 
multiple distance estimates can be potentially available for locating one 
UE in UMTS. In GSM this is not possible since the TA is available only 
for the unique serving cell. 

In the CI+RTT location method the unknown geographical coordinates of 
the UE whose position it is desired to determine are estimated by 
combining absolute distance measurements between the UE and the 
base stations in the active set. Each absolute distance measurement is 
calculated from each (RTT, RxTxTD) pair. 

In real applications it may happen that no <RTT, RxTxTD) pairs are 
available for location calculation. This may be due for instance to a 
measurement process failure or to a UE not supporting RxTxTD 
measurements. In such circumstances the CI+RTT location method -can 
still perform successfully a location calculation only on the basis of the 
knowledge of the CI of the serving cell {or of the multiple serving -cells if 
the UE is in soft handover). 

Even when all (or some) of the (RTT, RxTxTD) pairs are available, the 
location calculation based on the serving celKs) identity can be Used to 
improve performance of location algorithms that make use of RTT and 
RxTxTD measurements. 
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Thus the CI+RTT location method can use similar algorithms as the ones 
used by the SAI location method to estimate geographical coordinates and 
associated confidence region of a UE when the only location information . 
available is the identity of the serving cells. 

Having explained two location methods within which the invention can operate, 
embodiments of the invention will now be described, firstly with reference to 
figure 2. As mentioned previously, the embodiments of the invention involve 
the two main steps of : 

1 . estimating the location of the UE in terms of x - y coordinates and 

9 

2. calculating a confidence region for this location estimate. 

* 

A confidence region is a geographical region where the exact UE location is 
estimated to be with a given probability, referred to as the -confidence 
coefficient 0 < £ £ 1. 

Embodiments of the invention are implemented using the following location 
calculation methods (step 1 ) : 

a) Analytical Single-Cell Location Method 

b) Multi-Cell Location Method . 

c) Approximated Multi-Cell Location Method 

These location calculation methods are implemented by location calculation 
algorithms. Two classes of location calculation algorithms are used : 

• Location Estimate Calculation Algorithms to implement step 1 

• Confidence Region Calculation Algorithms*) implement step 2 

* 

Figure 2 is a flow chart showing the logical structure of embodiments of the 
invention. At the top of figure 2 is a high level location procedure 10 which 
controls all operations. The location procedure 10 calls lower level Location 
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Estimate Calculation Algorithms (LECAs) {step 1 ) and Confidence Region 
Calculation Algorithms (CRCAs) <step 2). LECAs calculate a location 
estimate by implementing one of the three available location methods a) to c) 
described below as three embodiments of the invention. 

Figure 2 shows that location estimate and confidence region calculation 
algorithms are kept logically separated. Thus below the location procedure 
10 are shown three location estimate procedures, an. "Analytical Single-Cell" 
procedure 12, a "Multi-Cell" procedure 14 and an "Approximated MuttnCell" 
procedure 16. Below the three location estimate procedures 12, 14, 16 are 
shown three confidence region procedures 18, 20, 22, corresponding 
respectively to the three location estimate procedures. An arrow from each 
location estimate procedure to its respective confidence region procedure 
shows that a confidence region is calculated only after a location estimate is 
calculated successfully. The possible ways of calculating a confidence 
region will be described below. 

When it is desired to obtain the location of a UE in the network 1, the 
location procedure 10 calls one of the three location estimate procedures 
12, 14, 16. The choice of location procedure is influenced by a user- 
defined variable TreferredLocationMethod". This variable indicates a 
preference towards a Location Calculation Method that should be used first 
for calculation of the location of the UE. The variable is set according -to 
the type of base stations in the SA and the concentration of them, with the 
aim of using the type of calculation that is most likely to succeed first. If 
the chosen location calculation method, fails for any reason, as shown in 
figure 2, control reverts to the location procedure 10, which decides 
whether to call a new LECA or to terminate the procedure with a failure. 
The location calculation could fail for a number of reasons, for example 
because for certain SA configurations the calculation according to a 
specific LECA might require an amount of resources (memory, 
computational capacity, etc.) that, at that specific moment, are not 
available in the system where the LECA is implemented. 
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When a LECA succeeds the CRCA calculating the confidence region 
according to the same Location Method implemented by the LECA that 
succeeded last is called. If the CRCA also succeeds the location 
procedure terminates successfully. On the other hand, if the CRCA fails, 
as shown in figure 2, control returns to the location procedure 10, which 
can then decide whether to try a different CRCA. 

The combination of location estimate and confidence region parameters is 
referred to as "shape". The shape definitions supported by the location 
calculation algorithms described above are: 

(i) Point Shape (i.e. including only the location estimate) 

(ii) Point with Uncertainty Ellipse Shape (where the 
confidence region is an ellipse) 

■ 

(iii) Point with Uncertainty Polygon Shape {where the confidence 
region is a polygon) 

(iv) Point with Uncertainty Arc Shape {where the confidence 
region is a circular arc) 

Diagrammatic representations of these four shape types are ^shown in 
figures 3a-3d. In the figures (Xest»ye$t)corresponds to (x,y) and x^m, yongm 

corresponds to xo.yo. 

Figure 3a shows the point shape. This has only one feature : 
• Coordinates of the location estimate x and y 

Figures 3b shows the point with uncertainty ellipse shape. This has the 
following features : 



WO 2004/028170 PCT/IB2003/004331 

12 

ft 

• Coordinates of an origin <at the location estimate) <x,y) 

• Semi-major axis R1 and semi-minor axis R2 

• Orientation angle a 

• Confidence coefficient value f 



Figure 3c shows the point with uncertainty polygon shape. This has the 
following features : 



• Coordinates of a location estimate (x.y) . 

• Number of vertices N {in the figure N = 7 ) 

• Coordinates of the vertices (xi,yi),. . . ,(Xn,Xn) 

Figure 3d shows the Point with Uncertainty Arc Shape. This has the 
fallowing features : 



Coordinates of a location estimate (x,y) 

Coordinates of an origin Xo and y 0 
Inner radius R1 and uncertainty radius R2 
Offset (orientation) angle a and included angle fi 
Confidence coefficient value f 



The CRCA used to determine the confidence region depends both on the 
location method used to determine the location estimate and on the shape 
type used to represent the resuiis. The shape type is influenced by a user- 
defined variable : "PreferredShapeType". This variable indicates a 
preference towards a shape type that should be used to represent the 
results. This variable -can be set appropriately to best show the-corifidence 
region in dependence on the characteristics of the SA. 

ft 

Finally, the choice of LECA and CRCA is influenced by a third user-defined 
variable "MethodsAndShapesAllowed". This variable defines a list of 
location method and shape type combinations that the location procedure 
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is allowed to use. Certain combinations may be incompatible, as some of 
the examples below will explain. 

A logical description of the location procedure is indicated in the following : 

Stage 1 : Create the list of Location Estimate Calculation Algorithms LECAList 

The list LECAList includes the Location Estimate Calculation Algorithms 
(LECAs) to be tried to determine a location estimate. The list is created by 
taking into account the following constraints: 

1 . Give highest priority to the LECA implementing the 
PreferredLocationMethod or, if only one cell must be considered, to those 
implementing the "Analytical Single-Ceir method. 

2. Add to the list only the LECAs implementing Location Methods that are 
allowed (i.e., those appearing at least once in the list 
MethodsAndShapesAllowed) 

Loop 1 :Loop over the list LECAList, trying only once each LECA in the list. 

As soon as one LECA in the list LECAList succeeds in determining a location 
estimate, go to Stage 2. If none of the LECAs in the list LECAList succeeds in 
determining a location estimate or if all LECAs in the list LECAList have 

» 

already been tried, go to Stage 4. 

Stage 2: Create the list of Confidence Region Calculation Algorithms 
CRCAList 

The list CRCAList includes the Confidence Region Calculation Algorithms 
(CRCAs) to be tried to determine a confidence region. The list is created by 
taking into account the following constraints: 

1 . Give highest priority to the CRCA that determines a confidence region with 
the shape specified by PreferredShapeType. 
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2. Add to the list only the CRCAs that are allowed (i.e.. those delivering a 
confidence region that, according to the list MethodsAndShapesAllowed, 
the Location Procedure is allowed to use when the LECA is the one that 
succeeded during Loop 1) 

Loop 2 : Loop over the list CRCAList, trying only once each CRCA in the list 

As soon as one of the CRCA in the list LECAList succeeds in determining the 
confidence region go to Stage 3. 

If none of the CRCAs in the list CRCAList succeeds, return to Loop 1 and try 
the next LECA in the list LECAList. 

Stage 3: Location Procedure Terminated with a Success 

The location result is represented in the shape type obtained by combining 
the latest location estimate determined inside Loop 1 and the corresponding 
confidence region determined inside Loop 2. 

Stage 4:. Location Procedure Terminated with a Failure 
A failure occurred because all LECAs in the list LECAList were tried and 
either none of them succeeded, or although some <or all) of them succeeded, 
none of the corresponding CRCAs in the list CRCAList succeeded in 
determining a confidence region. 

Having explained the logical process steps involved in embodiments of the 
invention, the LECAs and the CRCAs will now be described. 

LECAs 

The first consideration is modelling the geographical extension of each 
network cell of interest. The model is created by means of the following 
radio network parameters: 

• x - y coordinates for the Base Transceiver Station XBTS) antenna, 
(xs,ys) 
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• Bearing of the BTS antenna measured counterclockwise from x-axis 
in radians, <ps 

• Half Power Beam Width (HPBW) of the BTS antenna measured in 

~ - ■ * • 

radians, Aq>s 

• Maximum front radius of the cell, R?. This parameter specifies the 
maximum radius of the geographical region illuminated by the main 
radiation lobe of the BTS antenna, where the cell is serving. 

• Maximum back radius of the cell, Re. This parameter specifies the 
maximum radius of the geographical region illuminated by the back 
radiation lobe of the BTS antenna, where the cell is serving. 

« 

Figure 4 shows a representation of the cell border in the x-y plane. Based 
on the five radio network parameters listed above, each cell of interest is 
mathematically modeled in the (x,y) plane as follows : 

cit «\ - / " Rf ; 0 - W**) " M - A * 5 iii 

where d( x ,y) = J(z s -z) 2 + (y$-y) 2 d(x,y) and vix>y) is such that tan 

(pfr,y)= ys ~ y (d and if/ so defined are the radial and angular coordinates 
zs -z 

of a polar reference system centered at the BTS site). 

The border of the region covered by the Afceiis cells, S, is finally modeled by 
the algorithms as the union of Afceiis cell borders {Ci.... CNceiisJ-each of 
them modeled as in equation (1): 



fad 



<2) 



WO 2004/028170 PCT/IB2003/004331 

16 



The principles of the LECAs will now be discussed. The algorithms calculate 
the location estimate coordinates at the mass center of the geographical 
area covered by the Nce/fe cells of interest. According to this principle, the 
location estimate coordinates ( i.y ) are calculated using the following 

definitions: 

J J i*{x,y)xdxdv ^ J JT p(x,y)ydxdy 
J Ja(x,y)dxdy ' J j^x^dsdy ^ 

» 

where S is the border of the geographical region covered by the ficeiis cells 

> 

of interest and //(x,y) is the density of the geographical region enclosed by 
S. By density is meant the density of users, measured in units such as 

■ 

■ 

number of users per square km. 

In the context of the SAI location method S represents the Service Area 
(SA) border and p(x,y) the Service Area density. In embodiments of the 
invention, two alternatives for //<x,y) are available, according to the value 
of an algorithm parameter "-ConsiderOverlapping". This parameter can be 
set to TRUE or FALSE. The two alternatives are described in the following 
and exemplified in figure s. 

Alternative 1 : ConsiderOverlapping= FALSE 

The first alternative is to assume a uniform density over the whole SA. In 
this case (figure 5a), //(x,y) is defined as follows : 

* 

WW — i g elsewhere (4) 



where fj 0 is the constant density. Under this assumption the location 
estimate coordinates are expressed as 
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<5) 

where M(S) = H s dxdy is the area of the region confined by the Service Area 
border S. 



Alternative 2: ConsiderOverlapping= TRUE 

According to the second alternative, when ConsiderOverlapping= TRUE, 
p(x,y) is set to a value proportional to the number of cells covering the 
location of coordinates (x,y). In this case (figure 5b) the density is ^qual to 
the sum of Nc^s uniform densities of the cells in the Service Area: 



<6) 



where the density of the /-th cell /Ji{x,y) is equal to the constant fj\ Q over the 
cell area: 



W>y)-\Q elsewhere 

(7) 

Under this assumption the UE location estimate ^equation <3)) results as 

follows: 

\ 

Nam, * . Ncau , , 

L, fro I xdxdy 52 Pa ytedy 

i=l J J Ci j_i J Jd 
x — — — • t) — — 

"Cell. ' v Notllt 



i=l 



<8) 
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where, M(C i ) = jl Cj M{x,y)dxdy is the area of the /-th cell. 

From the comparison of this result with equation <5) it is evident that the 
location estimate (equation (8)) is the weighted average of the Ncells mass 
centers of the cells of interest; where the 

weight of the i-th mass center is ^oca» fj^M{Ci) = ' * " CeUs 



It should be remarked here that in examples described below, // 0 and //i 0 
are both set to 1 to simplify the calculations for illustrative purposes. 



Confidence Region Calculation 

As explained previously, the location estimate coordinates and the 
confidence region parameters resulting from the calculation performed by a 
location calculation algorithm are combined and represented with any of 
the four shape formats represented in figure 3. 

in practice all location methods except method a) "Analytical Single-Cell 0 
return a confidence region with the shape of a circle when they are 
requested to provide either a "Point with Uncertainty Arc" or a "Point with 
Uncertainty Ellipse" shape type. The circle is. centered at the Jocation 
estimate coordinates and has a radius Rcr, in general dependent on the 
confidence coefficient o <4 <> i . For this reason, in the following examples 
no difference is made, unless explicitly stated, between Arc and Ellipse 
confidence region calculations in the case of "Multi-Cell" and 
"Approximated Multi-Cell" location methods. 



It can further be understood that a circular confidence region is 
represented using the "Point with Uncertainty Arc" shape format by setting 
the inner radius Ri to zero, the uncertainty radius R2 to Rcr and the 
included angle p to 2k <the offset angle a is meaningless). See figure 3 
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for a graphical representation of these variables. The same circular 
confidence region is represented using the "Point with Uncertainty Ellipse 0 
shape format by setting both the semi-major axis /?i and the -semi-minor 
axis Rz to the radius Ron {the orientation angle a is meaningless). See 
figure 4 for a graphical representation of these variables. 

Examples 

First Embodiment- Analytical Single-Cell Algorithms 

The "Analytical Single-Cell" location method can be applied when only 
one cell is to be considered. In the following the algorithms implementing 
the location estimate and confidence region calculation are presented. 

Step 1 : Location Estimate Calculation 

The location estimate is calculated at the mass center of the unique cell 
of interest (Nceiis = 1). The coordinates of the mass center of the cell -can 
be calculated by evaluating analytically the integrals in equation (3). It-can 
be shown that by assuming a uniform density the following expression for 
the location estimate results: 



Step 2 : Confidence Region Calculation 
First variation : Ellipse 

* * 

The elliptical confidence region has the shape of a circle, the radius of 
which, Rcr, is calculated by scaling the distance between the location 
estimate and the furthest point on the cell borders (Rmax) by a factor 
equal to the square root of the confidence coefficient o < £ s l : 




<9) 
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Rcr = fa 



MAX 



tio) 



In this way the area of the circular confidence region having its centre at 
the location estimate and a radius equal to R C r has a total area equal to a 
fraction f of the circle centered at the same point enclosing the whole cell 
(i.e., corresponding to £ = 1). 

The maximum distance between the location estimate and the-cell 
borders is defined as follows : 

where R e ,b is the distance between the location estimate and the furthest 
point on the back cell region, and R £ , F is the distance between the 
location estimate and furthest point on the front cell region. These two 
distances can be calculated geometrically, as represented in figure6. Rb,b 
is defined as 



Re,b = V(* - x b¥ +(y~ y*f 



<12) 



where 




(13) 



R EF is defined as 



R B<P = <J(R F an&<j>s) 2 + (Rb + cos - Re#) 



04) 



Second Variation : Arc 
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The parameters of an Arc shaped confidence region can be calculated by 
means of trigonometric formulas. 

Third Variation : Polygon 

Given the maximum number of confidence region vertices, N, the polygonal 
confidence region is simply determined by generating N equally spaced 
pixels along the perimeter of the cell. The pixels correspond to the 
confidence region vertices. Since the polygonal confidence region has no 
confidence coefficient associated with it (see shape type definitions 
described with reference to figure 3) there is no need to take into account 
the confidence coefficient in this calculation. 

Second Embodiment : Multi-Cell Algorithms 

In the multi-cell location method a location estimate is determined 
numerically at the mass center of a uniform rectangular grid covering the 
geographical region enclosed by the multiple cell borders. The method 
can be applied when one or more cells is to be considered. 

The rectangular grid, which is used also to determine the confidence 
region, is obtained by sampling uniformly in the x and y directions the 
area covered by the cells of interest, using constant step sizes ax and Ay, 
respectively. An example of a grid created is shown in figure 7a, and is 
defined by Np sets of three values: 

♦ 

<15) 

where 

• (Xp,y p ) are the central x - y coordinates^ the p-th rectangular 
pixel having area A % * Ay; 

• w p is the number of cells that overlap in the area represented 
by the p-th pixel. 
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Step 1 : Location Estimate Calculation 

The location estimate is determined through a discretised version of 
equation (3): 



N P Np 



(16) 



where 



(17) 



is the normalised density associated with the p-th pixel, assumed 
constant over the pixel area. 

Step 2 : Confidence Region Calculation 

First Variation : Circle 

Given the location estimate coordinates (x.y) the following set of 
distances between the location estimate and the grid pixels <Xp,y p } £p, 
calculated : 

The p-th distance is then weighted by the normalised density of the 

■ 

corresponding pixel, q p . and the resulting set of distances Np and 
corresponding weights is arranged in a discrete distribution defined by the 
following Nd z Np sets of two values: 



i 



■ 
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♦ 

vhere d x < <fe < . . . < awb are the re-arranged Afo distances and {Pj}"i are 
the probabilities of the N D distances such that = 1. 

Subsequently, the radius of the circular confidence region Rc/?is calculated 
as the (100£)-th percentile of the distance distribution <0 < 4 s 1): 

JO 

i?CR — djo with jo such that J^Pi ~ £ 

*=i (20) 

Figure 7a shows the circular confidence region obtained using the 
algorithm described above on a three-cell sample using a confidence 
coefficient 4 = 0.8. The same figure shows the pixel coordinates and 
location estimate as well. Figure 7b shows the discrete distance distribution 
(equation (19)) and the corresponding 80-th percentile. 

Second Variation : Polygon 

Given the maximum number of confidence region vertices, N, the polygonal 
confidence region is obtained from the grid of Np pixels defined in equation 
(15) by properly selecting N outermost pixels that best approximate the grid 
borders. Since the polygonal confidence region has no confidence 
coefficient associated with it <see shape type definitions described with 
reference to figure 3) there is no need to take into account the confidence 
coefficient in this calculation. 

Figure 8 shows the polygonal confidence region resulting from the same 
sample cells used to produce the results shown in figure 7. 

* « 
Effect of Service Area Density Assumption 
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In one implementation of the multi-cell location estimate calculation 
algorithms the density //(Xp.Yp) used in equation (17) is defined according to 
the value of the parameter ConsiderOveriapping. When 
ConsiderOverlapping is FALSE the number of cells overlapping on each 
pixel area is not considered when -calculating the location estimate, but 
when ConsiderOveriapping is TRUE, the overlapping is considered. These 
assumptions result in the following definition forp(Xp,y p ): 

J 1 ConsiderOveriapping = FALSE 
WWW ~\w p ConsiderOveriapping = TRUE (21 ) 

In the implementation when ConsiderOveriapping is TRUE the location 
estimate is not calculated with equation (16) but using explicitly equation 
(8), which is only approximated by equation {1 6). 



Third Embodiment : Approximated Multi-Cell A lgorithms 
In this embodiment a location estimate is determined numerically at the 
mass centre of a polygon P enclosing the borders of the multiple cells. The 
method can be applied when one or more cells must be considered. 

The polygon P is obtained as a polygon of Ny vertices enclosing pixels 
equally spaced along the borders of the ceiis of interest For location 
calculation purposes, the polygon P is used as an approximation of the 
border enclosed by the multiple cells (P approximates the bonier S 
defined in equation (1)). Figure 9 shows an example of how a polygon 
enclosing three cells is obtained from pixels placed on the borders of the 
cells. 

Step 1 : Location Estimation Calculation 

The location estimate coordinates <x,y) are calculated at the mass center 
of the polygon P approximating the cells' border. Analytical formulas can 
be used for this purpose. 
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First Variation : Circle 

The radius of the circular confidence region Re* is calculated by scaling the 
distance R M ax between the location estimate and the furthest vertex of the 
polygon P by a factor equal to the square root of the confidence coefficient 
0< 5 <; 1: 

RCR — ifZRhfAX (22) 

In this way the area of the circular confidence region having its centre at 
the location estimate and a radius equal to Rcr has a total area equal to a 
fraction f of the circle centered at the same point enclosing the whole 
polygon P used to approximate the border of the multiple cells (i.e., 
corresponding to £ = 1). Figure 10a shows the circular confidence region 
resulting from the above calculation applied to the sample cells of figure 9 
assuming 4 = 0.8. 

Second Variation : Polygon 

Given the maximum number of the confidence region vertices, N s if N < N V 
the polygonal confidence region is obtained by properly selecting N vertices 
among the N v vertices of the polygon P. If A/ s N v , the polygonal confidence 
region coincides with P itself. Since the polygonal confidence region has no 
confidence coefficient associated -(see shape type definitions as described 
with reference to figure 3) there is no need to take into account the 
confidence coefficient in this calculation. Figure 10b shows a polygonal 
confidence region resulting from the sample cells of figure 9 ('m the figure N. 
= 15). 

* 

Effect of Service Area Density Assumption 
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In one implementation of the third embodiment the polygon P is a convex 
polygon. The implementation can be effected with the parameter 
ConsiderOverlapping. If ConsiderOverlapping is set to 
FALSE the calculations of this Third embodiment as 
described above are performed. If, on the other hand, 
ConsiderOverlapping is set to TRUE the method is in practice the 
"Multi-Cell" location method (the second embodiment) with the difference 
that the grid (equation (15)) instead of being a uniform rectangular grid is a 
non rectangular grid obtained by subdividing each cell in the Service Area 
into Np circular arcs having approximately the same area, and by 
assigning : 

. the x - y coordinates of the p-th pixel <Xp,y p ) at the center of the p-th 
circular arc; 

• the weight w p equal to the area of the p-th circular arc. 

An example of non rectangular grid created with this method and covering the 
area of one sample cell is represented in figure 1 1 . 

The applicant hereby discloses in isolation each individual feature described 
herein and any combination of two or more such features, to the extent that 
such features or combinations are capable of being carried out based on the 
present specification as a whole in the light of the common general knowledge 
of a person skilled in the art, irrespective of whether such features or 
combinations of features solve any problems disclosed herein, and without 
limitation to the scope of the claims. The applicant indicates that aspects of 
the present invention may consist of any such individual feature or 
combination of features. In view of the foregoing description it will be evident 
to a person skilled in the art that various modifications may be made within 
the scope of the invention. 
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